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Abstract 

We compute the Bk parameter and the mass difference AM^s of the K°-K° system by 
means of the chiral quark model. The chiral coefficients of the relevant AS* = 2 and 
AS = 1 chiral lagrangians are computed via quark-loop integration. We include the 
relevant effects of one-loop corrections in chiral perturbation theory. The final result is 
very sensitive to non-factorizable corrections of 0(asN) coming from gluon condensation. 
The size of the gluon condensate is determined by fitting the experimental value of the 
amplitude K + — ► 7r + 7r°. By varying all the relevant parameters we obtain 

B K = 0.87 ±0.33 . 

We evaluate within the model the long-distance contributions to AMis induced by the 
double insertion of the AS = 1 chiral lagrangian and study the interplay between short- 
and long-distance amplitudes. By varying all parameters we obtain 

AMl h s /AM e L f = 0.76 t°oM ■ 

Finally, we investigate the phenomenological constraints on the Kobayashi-Maskawa pa- 
rameter Im At entering the determination of e' je. 
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The mixing in the K°-K° system is determined by the weak effective hamiltonian 
through the mass matrix (K°\TLw\K°) . Its computation requires both long- and 
short-distance physics. For this reason, the CP violating quantity e and the mass 
difference AM^s, which are respectively related to the imaginary and the real part 
of the mass matrix, are factorized in the renormalization-group invariant parameter 
E>k, which takes into account long-distance effects, and the Wilson coefficients, 
which account for the short distance ones. AM^g receives in addition a long-distance 
correction from the double insertion of the AS = 1 hamiltonian. 

In this work we apply the chiral quark model []J (xQM) and chiral perturbation 
techniques to estimate the long-distance part of the mass matrix. Such an analysis 
follows our recent studies of the AI =1/2 selection rule || and e'/e || within the 
xQM approach to kaon physics. 

The long-distance contributions thus computed are eventually matched to an up- 
to-date next-to-leading order (NLO) determination of the short-distance coefficients 
and the results compared with the experimental values. 



1 Introduction 



The effective AS = 2 quark lagrangian at scales fi < m c is given by 
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(1.1^ 



where Gf is the Fermi constant, Mw is the W boson mass, Xi = mf/M^, and \x 
is the renormalization scale. The parameters Aj = VjdV* s represent the relevant 
combinations of Kobayashy-Maskawa (KM) matrix elements (j = u,c,t). Finally 
we denote by Qs2 the AS = 2 local four quark operator 

Qs2 = {s L ^d L ){s Ll ^d L ). (1.2) 

The integration of the electroweak loops leads to the Inami-Lim functions 
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which depend on the masses of the charm and top quarks and describe the AS = 2 
transition amplitude in the absence of strong interactions. 

The short-distance QCD corrections are encoded in the coefficients rji, rj2 and 773 
with a common scale-dependent factor b(fi) factorized out. They are functions of 
the heavy quarks masses and of the scale parameter Aqcd- These QCD corrections 
are available to NLO || |], [7[] in the strong and electromagnetic couplings. 

The scale-dependent common factor of the short-distance corrections is given by 

%) = [«,(/i)]- 2/9 (l-J 3 ^) , (1.5) 

where J3 depends on the 75-scheme used in the regularization. The naive dimensional 
regularization (NDR) gives 

7 ndr _ (ir) 
J3 - 162 , U-bJ 

while in the 't Hooft-Veltman (HV) scheme one finds 

7 hv _ (1 7^ 

J3 ~ 162' [L7) 

For the running QCD coupling we take the average over recent LEP and SLC 

determinations ||, 

a s (m z ) = 0.119 ±0.006 (1.8) 
which corresponds to 

Aqq D = 350 ± 100 MeV . (1.9) 
The scale-dependent Bk(/J<) parameter is defined by the matrix element 
{K»\Q S2 (//) \K°) = p 2 K m 2 K B K (f,) , (1.10) 

where fx and tuk are the kaon decay constant and mass, respectively (see table 2 
for their numerical values). 

The value of S^(/i) measures the deviation of the matrix element from the 
vacuum saturation approximation used in the original work of Gaillard and Lee , 
namely Bk{^) = 1. The physically relevant parameter is Bk, which is defined by 
the relation: 

B K = B K (n)b(ji) . (1.11) 
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B K 


Method 


Reference 


3/4 
0.37 


Leading 1/N C 
Lowest-Order Chiral Perturbation Theory 


1 


11 

12 


1 


0.70 ±0.10 

0.4 ± 0.2 
0.42 ±0.06 
0.60^0.80 


Next-to-Leading 1/N C Estimate 
Next-to-Leading 1/N C Estimate, 0(p 2 ) 
0(p 4 ) Chiral Perturbation Theory 
NJL model with spin-1 interactions 


1 


13 
14 
15 
16 


1 
1 

1 


0.39 ±0.10 
0.5 ±0.1 ±0.2 
0.55 ±0.25 
0.58 ±0.22 


QCD-Hadronic Duality 
QCD Sum Rules (3-Point Functions) 
QCD Sum Rules (3-Point Functions) 
Laplace Sum Rule 


1 
1 


7 

19 

20 
21 


1 
1 


0.90 ±0.03 ±0.14 


Lattice 




22 





Table 1: Values of B K obtained in different approaches. 



This quantity should be in principle renormalization scale independent. As we 
include the perturbative NLO determination of the Wilson coefficient, we shall also 
discuss the 75-scheme dependence of our result. 

A useful up-to-date summary of various determinations of this parameter is given 
in Table 1 which updates that of ref. [|l(| . 

We have followed the approach described in ref. |23| in which the weak chiral 
lagrangian is considered as the effective theory of the xQM |1| . In the present case, 
it is the bosonization of the operator Q52 and the determination of the coefficient 
of the corresponding AS = 2 chiral lagrangian that is made possible by the xQM. 

In the determination of Bx(fi) to 0(a s N c ) enters the contribution of the gluon 
condensate. The final estimate is very sensitive to the value of such an input pa- 
rameter. In order to restrict the range of allowed values, we impose the additional 
constraint of taking for the gluon condensate the value that gives the best fit of the 
experimental amplitude K + — > 7r + 7r°, which is related at the leading order in chiral 
perturbation theory to that of K° — > K°. Such a procedure is consistent with that 
followed in ref. || where we reproduced the AI = 1/2 rule by a similar choice of 
input parameters. 
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A "long-distance" scale dependence is introduced by the one-loop chiral cor- 
rections to the hadronic matrix elements. According to our approach, this scale 
dependence should match that in the Wilson coefficients and provide a scale inde- 
pendent value of Bk within the uncertainties quoted for higher order corrections. 
In the case of Bk, we find that there cannot be matching at this order insofar as 
both the Wilson coefficient and the chiral corrections renormalize the parameter in 
the same direction. The scale dependence remains however below 20% and the final 
estimate is thus still reliable. 

Our approach is in principle sensitive to the scheme used to treat 75 matrices 



in a generic space-time dimension. The NDR prescription of ref. [gj, |25] preserves 
the chiral properties of the operator Qs2 by means of a convenient normalization 
of the evanescent operators. As discussed in ref. P^| , the consistency with such a 
prescription makes the matrix elements of Q52 the same in the two schemes. As a 
consequence, the remnant scheme dependence of the final result is that present in 
the short-distance factor b(/i). 

There are two important parameters related to the K°-K° mixing: the CP 
violating quantity e which is proportional to the imaginary part of the mass matrix 
and the mass difference AM LS = m L — m S - The observed value for these quantities 
are [^6| : 

|e| = (2.266 ± 0.023) x 10~ 3 (1.12) 

and 

AM LS = (3.510 ± 0.018) x 10~ 15 GeV . (1.13) 

Knowing e, we can determine ImA t , as discussed in section 6. As a by-product of 
the computation one also obtains an estimate for the width difference ATls, the 
experimental value of which is 

AT LS = -(7.374 ± 0.010) x 10~ 15 GeV . (1.14) 

However, a consistent determination of this quantity requires one extra order in 
perturbation theory, as we shall discuss below. 

From the theoretical point of view, the K° — K° mass matrix can be written, 
using CPT invariance, as 

1 / (K°\H W \K°) (K°\H W \K°) 



M 



2m K \ (K°\H W \K°) (K°\H W \K°) 
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In the presence of CP violation (e ^ 0) M 12 and T 12 are complex numbers. The 
diagonalization of the mass matrix ( |1.15|) leads to the physical states: 



K L 
K s 



V20+FF) 
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K°(l + e)-K°(l-e 
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For a tiny CP violation, their associated mass and width differences are given by: 
1 
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In order to estimate these two parameters we need to evaluate in addition to the 
quark box-diagram contribution, coming from the AS = 2 effective weak lagrangian 
given in (1J-), the long-distance contribution coming from the double insertion of 
the AS* = 1 weak chiral lagrangian. In the latter case, the mixing between K° and 
K° can proceed, up to the one-loop level, via one- and two-particle intermediate 
states 



;i.i9) 
;i.2o) 



Within the xQM approach the AS" = 1 weak chiral lagrangian can be systemat- 
ically derived at a given order in momentum expansion starting from the effective 
quark lagrangian P7| : 



AS=1 



;i.2i) 



where Qi are local four-quark operators obtained by integrating out in the standard 
model the vector bosons and the heavy quarks t, b and c. A convenient and by now 
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standard basis includes the following ten quark operators: 

Qi = (s a up)y_ A (u/3(i a ) v _ A , 
Q 2 = (sm) v _ a (ud) y _ A , 
Q 3 ,5 = O0 v _ A £ ? (gg) VTA , 

Q 7j9 = | (sgT) v _ a Eg e g (gg) v±A , 

Qs.io = | (s Q <i j a) v _ A E g eq(g / 3g Q )v±A , 

where a, [3 denote color indices (a, (3 — 1, . . . , iV c ) and e q are quark charges. Color 
indices for the color singlet operators are omitted. (V ± A) refer to 7^(1 ± 75). We 
recall that Q\ i stand for the ^-induced current-current operators, £^3-6 for the 
QCD penguin operators and Q7-10 for the electroweak penguin (and box) ones. 

The functions and yi(fi) are the Wilson coefficients and Vij the KM matrix 
elements; r = -V td V*/V ud V* s . 

In a previous work |23| we have computed the chiral coefficients for the complete 
0(p 2 ) AS = 1 chiral lagrangian. We will make use of those results to evaluate the 
long-distance contributions in eqs. ( |1 . 19| ) — ( |1 ■ 2Q[ ) . 



2 A MODEL INDEPENDENT ESTIMATE OF B K 

The AS" = 2 matrix element can be related via chiral symmetry to that of the 
AS = 1 and AI = 3/2 amplitude A(K+ -> 7r+7r°)[IJ. Neglecting the SU(3) 



breaking effects related to the chiral loop corrections to the matrix element, the 
electromagnetic contributions and the 77 — 77 mixing, we obtain the relation 

\fWA = p. J- ,f« , - M A (K+ *V) . (21) 
3 G F V* s V ud m z K - mi z 1 (//) + z 2 (n) 

In the previous equation V us and V ud are two matrix elements of the KM mixing 
matrix, is the AS" = 2 Wilson coefficient given by (|1.5| ), while zi(fi) and z 2 (/x) 
are the real parts of the Wilson coefficients for the two AS" = 1 operators Qi and 
Q2 which dominate the K + — > tt+tt transition. 

By inputting the experimental value A(K + — > 7t + tt°) = 1.84 x 1CT 8 GeV and the 
NLO results for the Wilson coefficients (the ratio b(fi) / (zi(fj,) + z 2 {^)) is to a large 
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extent \i and 75-scheme independent) we find the model "independent" estimate 
B K = 0.40 . (2.2) 



This number updates the value Bx = 0.33 given in ref. [E^ . 

On the other hand, having a model that reproduces the experimental result, in 
order to apply correctly eq. ( |2.1|) we must subtract in A(K + — > 7r + 7r°) all the chiral 
symmetry breaking corrections due to chiral loops, electroweak penguins and tx — rj 
mixing ||. In this way we obtain in the xQM approach, on the basis of chiral 
symmetry arguments alone, the following 0(p 2 ) prediction: 



Bk = (ji) Ijf 
4 Ik 



1 + ^(1-^) 



(2.3) 



iv c \uc s ^^l n/ 

(GG) = — TT^TT ' ( 2A ) 



In the previous formula we have denoted by <5(gg) the non-perturbative gluonic 
corrections which arise in the xQM approach, 

N c ( a s GG/n } 

lQTT 2 f 4 

where (a s GG/7r) is the gluon condensate and iV c is the number of colors. We will 
come back to these corrections in the next section. 

In considering eq. ( |2.3|) it is important to remember that the factor f n comes 
from the soft pion theorem, while / is the chiral lagrangian parameter appearing in 
the calculation of the amplitude A(K + — > 7r + 7r°). At the tree level / = f n . The 
spurious /i dependence present in eq. ( |2.3| ) should be canceled by that of the hadronic 
matrix elements, which is absent at the lowest order in the chiral expansion. 

If we choose for the gluon condensate the value (a s GG/n) = (360 MeV) 4 (which 
gives the best fit of A(K + — > 7r + 7r )), we obtain at /i = 0.8 GeV 

B K ~ 0.33 . (2.5) 

This value includes the non-factorizable effects of gluon condensate corrections, 
which play a crucial role in the fit of the AI = 3/2 amplitude in K — > txtx de- 
cays. 

The value in eq. ( |2.5| ) represents the starting point of our analysis, to which we 
will add the effect of chiral loop contributions to the AS = 2 matrix element. 
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3 Computing B k 



In this section we will extend the techniques that we have developed for AS = 1 
weak processes in ref. p3[] , by using the xQM to construct the AS" = 2 weak chiral 
lagrangian. 

3.1 The leading chiral coefficient in the chiral quark model 

At the leading 0(p 2 ) order in the chiral expansion, the strong interaction between 



the SU(3) Goldstone bosons is described by the following effective lagrangian |f28 

4t 2 L g = y Tr (D^D^) + ^-B Tr (MY? + X-M 1 ") , (3.1) 
whereA^ is the mass matrix of the three light quarks (u,d and s) and E is defined 



as 



S = exp (y U(x) ) , U(x) = ^ a (x)/2 . (3.2) 

V / / o=1..8 

To the same order, the AS* = 2 weak chiral lagrangian is given by: 

4 2 L 2 = C(Qs2) Tr (XlYD^) Tr (\ 3 2 ED^) . (3.3) 



In eq. ( p.3|) indicates the covariant derivative with respect to any external field, 



while is a combination of the SU(3) Gell-Mann matrices which acts in the flavor 
space causing a transition from a <i-quark to an s-quark: (A^/fc = S 3 iS 2 k- 

C(Qs2) is the chiral coefficient, which we determine by comparison with the xQM 
calculation. Two configurations contribute to the determination of this coefficient 
at 0(N C ), as shown in Fig. 1. 

In both HV and NDR schemes we find (for convenience we do not write the 
overall AS* = 2 Wilson coefficient given in eq. ( |1 . 1| ) ) 

C{Q S2 ) = -£.(l + ±.). (3.4) 



4 V N c 

An important correction to eq. (|3.4| ) arises by considering the propagation of 
quarks in an external gluon field. The effects of non-perturbative gluonic corrections 



have been first studied in [|i4~ 



In the case of K° — > K° transition the relevant gluonic corrections are given 
by the diagrams of Fig. 2. By including gluonic condensate corrections, eq. 
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K° 



— <- 




K° 



--<- 



Figure 1: The two configurations relevant to the determination of the chiral coeffi- 
cient C(Qs2)- The crossed circles stand for the insertion of the currents of the local 
AS = 2 four-quark operator. 




K° K° 
«<--- -->■ 




K 1 



Figure 2: Same as Fig. 1 with the inclusion of the gluon corrections. 



becomes 
C(Qs2) 



f 4 r i 
J i + — (l 

Nr. 



\GG) 



(3.5) 



where 5(gg> is given by eq. (|2.4j). 

By using the definition given in eq. ([L.10|) and computing at leading order 
(K°\Qs2\K°), we obtain the following expression for B K (/i): 



Bk(p) = j 



Hgg) 



P 
f 2 



(3.6) 



K 



At this stage of the computation, Bk(h) does not exhibit yet an explicit dependence 
on jj,. In our approach the scale dependence arises from meson-loop corrections. 
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\ / 
K W' K 
-->-■—<-- 

a 



V 



— > 



K WK° — K° 
—>-»->--■-<— + 

Figure 3: One-loop chiral corrections to the kaon mass matrix. The black box and 
circle indicate the insertion of the weak AS = 2 and strong chiral hamiltonians 
respectively. The octet mesons K, it, rj are exchanged in the loop. 



If we take / = /„. in eq. Q3.6|) we recover eq. fl2.5|) , as it should be. 
Taking / = f K and 5{gg) — 0, eq. ( |3.(j| ) reproduces the result obtained in the 
1/N C approach. 



3.2 One-loop renormalization of the AS = 2 transition 

So far we have ignored chiral- loop corrections to the evaluation of Br. The 
introduction of these contributions gives a long-distance \i dependence to B^ip). 

In conventional chiral perturbation theory the scale dependence of meson loops 
renormalization is canceled by construction by the 0(p 4 ) counterterms in the chiral 
lagrangian. In our approach, on the contrary, the tree-level counterterms are // in- 
dependent and the scale dependence introduced in the hadronic matrix elements via 
the meson loops, evaluated in dimensional regularization with the standard minimal 
subtraction, is matched with the scale dependence of the Wilson coefficients. 

The diagrams relevant to the present case are depicted in Fig. 3. The diagram in 
Fig. 3(a) contains only a four- meson weak vertex of the AS = 2 chiral Lagrangian, 
while the diagram in Fig. 3(b) contains two vertices, one of which is a four- meson 
strong vertex and the other is a two-meson weak vertex. Another class of diagrams, 
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which are induced by wave function renormalization, is shown below 3(a,b). 
A direct calculation yields: 



(K°\Qs2\K°)\ tree = C(Q 



S2, 



-Am 2 K 



P 



(K°\Q S 2\K C 



(K°\Qs2\K° 



la 



C(Qs2) ji [3/4 (mj) + h (ml) + 7m 2 K I 2 (m 2 ) 



16 



+ 5m 2 K I 2 (ml) + -m 2 K I 2 (m 2 K ) 

h (m 2 K ) + 3m 2 K I 2 ( 



2 C(Qs2) -ji 
+ 3m K I 3 (m 2 K ) 



(K \Q S2 \K°)\ wf 



-4 C(Q 



m 



S2) 



A 



-h (ml) + \l, 



A' 



nr.. 



+ ~h (m 2 K ) 



where C{Qs2) is the chiral coefficient given by eq. fl3.5|), and 



h (m 2 ) 
h (mt) 
h (ml) 



(27T) 



(q 2 - ml) 
1 



(2n) 4 J (q 2 -m 2 ) 2 



d 4 q 



d 4 q 



1 2 ( , m i 



' lnfe 



- 1 



(2nfJ (q 2 ~m 2 ) 



d 4 q 



16tt 2 
f 



16n 



m A In 



2 ""1 



- 1 



(3.7) 

(3.8) 
(3.9) 

(3.10) 

(3.11) 
(3.12) 
(3.13) 



We also have to consider the meson decay constant renormalization, that is the 
one-loop determination of / in terms of fx' 



f = f K 1 + 



2/2 



-Mrnl) + -h{m K ) + -I 2 {m 2 ) 



(3.14) 



where the dots represent contributions of the 0(p 4 ) chiral lagrangian. This renormal- 
ization introduces chiral corrections which cancel some of the contributions coming 
from the meson loops in Fig. 3. Notice that the divergent integrals in eqs. ( J3.11p - 
( |3.13D are minimally subtracted, while the usual chiral expansion prescription keeps 
only the chiral logarithms. We will discuss the implications for the 0(p 4 ) countert- 
erms in a forthcoming publication |2^ . 
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The expression for the AS = 2 amplitude, comprehensive of chiral loops, wave 
function and kaon decay constant renormalizations can be written as 



(K°\Q S2 \K°) = (K°\Q S2 \K°)\ tree [l + 2 



f-l 



K 



K 



+ (K°\Qs2\K°) 



|la+16+w/ j 



(3.15) 



where the chiral lagrangian coefficient / in the term (K \Qs 2 \K )\ tree oc f 2 is re- 
placed by the physical decay constant fx (all other terms are independent on /). 
As it was noted by Bruno in ref. this replacement corresponds to resumming 
all orders of the chiral expansion in the factorizable component of the amplitude. 
An approach similar to the one we are adopting here has been followed in ref. 



| 13fl in the framework of a cut-off regularization of the chiral loops. It is important 
to stress that we have chosen to regularize the divergent integrals appearing in the 
meson loops by using dimensional regularization (as we have already done in [^. 
|[]). This choice is motivated by consistency with the short distance calculation of 
the Wilson coefficients, which is performed using the same regularization. 

In order to show the impact of chiral loops on the K°-K° amplitude, we find 
convenient to factorize the tree level contribution in terms of the relevant parameters, 
while giving the corresponding loop renormalization as a numerical coefficient with 
an explicit \x dependence. The values of the meson masses and other input variables 
are those given in Table 2. 

We thus find: 



(K°\Q S2 \K°)=m 2 K l 



K 



1 



Wc (i - s m 



1 + 0.728 + 0.372 In/ , (3.16) 



where /i is taken in units of GeV. From eq. (|3.16 ) we obtain the final result for 
Bk{h)i inclusive of the effects of meson loops, wave function and kaon decay constant 
renormalization: 



1 



i 



><GG> 



1 + 0.728 + 0.372 ln/r 



(3.17) 



The scale dependence of the hadronic matrix element interferes constructively 



with that of 6(/i) in eq. ( |1.5| ). Nevertheless, the overall scale dependence of Bk 
remains below 20% in the range between 0.8 and 1 GeV. 
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parameter 


value 


fir fir+ 


92.4 MeV 


fx = fic+ 


113 MeV 


m n = (m w + +m 7r o)/2 


138 MeV 


rriK = m K o 


498 MeV 




548 MeV 


a(4) 
iv QCD 


350 ± 100 MeV 



Table 2: Table of the numerical values used for the input parameters. 

4 Numerical Analysis 

We now have all the ingredients necessary to make a detailed analysis of the 
values of the parameter B^, where Bjc(^i) and are given by eq. ( [3.17]) and 
eq. ( p..5|) , respectively. The final result depends on the values of the gluon con- 
densate (a s GG /tt) entering in the determination of the gluon corrections to B k (ii) 
and of Aq^, d which determines the value of the QCD coupling constant as, and 
consequently of b(fi) . 

4.1 Input Parameters 

A relevant input parameter in our present analysis is the gluon condensate. We 
choose for this quantity the value that gives within a 30% error a fit of the A J = 3/2 
K + — > 7T + 7r° amplitude: 

{^-Gg\ = (360 ± 15 MeV) 4 . (4.1) 

Although this value fits leading order QCD sum rule determinations, the relation be- 
tween the chiral quark model gluon condensate and other estimates (e.g.: QCD sum 
rules) is rather unclear (only low frequency modes of the gluon fields are included 
in eq. (^4.1| )). We stress that our approach is in this respect purely phenomeno logi- 
cal: we consider the quark and gluon condensates as parameters of the model to be 
consistently determined by comparison with known observables. A redundant set 
of determinations provide the basis for a predictive framework. A word of caution 
concerning the renormalization prescription of the chiral lagrangian parameter / in 
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the amplitudes: in refs. |23|, |2|, |3|] we have included the one-loop renormalization of 
l// 3 in the K — > tttt tree level chiral amplitudes. From now on we include in the 
counting of powers of / also the / dependence of the chiral coefficient computed 
in the xQM. The numerical consequences of this change in prescription are that 
the best fit of the AI = 1/2 rule leads to a central value of the gluon condensate 
(a s GG/ir) = (360 MeV) 4 , slightly smaller than that obtained in 0, namely (372 
MeV) 4 , and to a central value of the quark condensate of (— 280MeV) 3 , slightly 
larger than (— 271MeV) 3 , quoted as the best fit in ref. 0. Since our present results 
depend very little on the quark condensate we keep it fixed at the value (qq) = — 
(280 MeV) 3 , while varying the gluon condensate in the range of eq. (fOj). 

The present analysis, as our previous ones, includes 0(p 2 ) chiral coefficients 
and the effects of chiral loops. A complete 0(p 4 ) calculation is under way [29 



Preliminary results show that the best fit of the AI = 1/2 rule is obtained for lower 
values of the quark and gluon condensates. 

Another input parameter which is important for the determination of Bk, is 
the QCD running coupling constant a s entering in the computation of the short 
distance factor 6(/u). In our numerical estimates we use for a s the range of eq.flLSp, 
corresponding to the values of Aq^ d given by (|1.9|). The values of this and other 
input parameters are listed in Table 2. 

4.2 Numerical results for Bk 

Our numerical estimate of the parameter Bk is summarized in Table 3, in which 
we have fixed (a s GG/ir) to the central value of eq. ( f4.1| ) and we have examined two 
extreme values of the matching scale /i in both schemes HV and NDR. The three 
parts of the table show the dependence on the QCD scale parameter Aq^d- 

From Table 3 we obtain the ranges 0.79 < B K < 1.0 in NDR and 0.69 < B K < 
0.97 in HV scheme. 

The quantities A^Bk and A^Bk measure the size of the 75— scheme and //— 
dependences respectively, 



A J5 B K 



A^Bk 



B k 


HV — B K 


NDR 




B k 


hv + Bk 


NDR 




B k 


(1 GeV) - 


-B A - (0.8 GeV) 



B K (1 GeV) + B K (0.8 GeV) 



(4.2) 
(4.3) 
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A OCD = 250 MeV 




/i = 0.8 GeV 


H = 1 GeV 




NDR 


HV 


NDR 


HV 


Mm) 


1.25 


1.19 


1.30 


1.24 


Bk 


0.88 


0.84 


1.01 


0.97 


A 15 B K 


5.2% 


4.2% 


A^B K 


14% - 15% 




4% - 5% 


Aqcd = 350 MeV 




fi = 0.8 GeV 


A» = 1 GeV 




NDR 


HV 


NDR 


HV 


b{p) 


1.17 


1.08 


1.23 


1.17 


B K 


0.83 


0.77 


0.96 


0.91 


A 75 B K 


8% 


5.7% 


A^B K 


15% - 17% 




5% - 7% 


d = 450 MeV 




fj, = 0.8 GeV 


/x = 1 GeV 




NDR 


HV 


NDR 


HV 


b{p) 


1.12 


0.97 


1.18 


1.09 





0.79 


0.69 


0.92 


0.85 


A 75 Bk 


14.1% 


7.9% 


A^B k 


15% - 21% 


A„6(/x) 


5% - 11% 



Table 3: Matching scale and 75 scheme dependence of Bk in the xQM with NLO 
Wilson coefficients, for various values of Aq^ d . We take for the gluon condensate 
the value (a s GG/ir) = (360 MeV) 4 , preferred by the fit of T(K+ -> tt+tt ). 
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I 

PQ 



0.8 



0.6 



0.4 




0.345 



0.35 



0.355 



0.36 
<GG> 



0.365 



0.37 



0.375 



Figure 4: The B K parameter is shown as a function of the gluon condensate for 
Aq^ d = 350 MeV and fi = 0.8 GeV. We denote by (GG) the quantity (a s GG/n) 1/4 
in units of GeV. The dark and grey lines represent the HV and NDR results respec- 
tively. 



The scale dependence A^B K is near to 20% in both schemes and it is mainly 
due to the effect of meson loops renormalization. As a matter of fact the final \i 
dependence is larger than the one originally present in the coefficient b(/i), which is 
less than 10%. Nevertheless the fact that the scale dependence is at most 20% makes 
us confident on the stability of our results, and allows us to choose \i = 0.8 GeV 
as the best compromise between the upper limit of validity of chiral perturbation 
theory, used to compute Bk((J>), and the lowest scale for perturbative calculations, 
needed to obtain the short distance coefficient b(fi). 

The scheme dependence of our result is entirely due to b(/j), since the hadronic 
matrix element does not exhibit any scheme dependence. At any rate A 15 Bk is 
below 10% for all values of \x in the given range. 

Finally, a few words on the dependence of our results on the value chosen for the 
gluon condensate. In Fig. 4 we show B^ as a function of the gluon condensate, for 
our preferred matching scale [i = 0.8 GeV, and Aqq D = 350 MeV. 

It appears that Bk is a sensitive decreasing function of the gluon condensate. 
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For our values of {a s GG/n) the term 5(gg) i n ec i- (|3.17|) is greater than 1, thus 



changing the sign of the 1/N C contribution and determining a reduction of the final 
result for Bk- By varying also the value of the gluon condensate in the range of 
eq. (|4.1|) we obtain the overall range 



0.54 <B K < 1.2 . (4.4) 

which represents our conservative prediction for Bk- The central value Bk = 0.87, 
quoted in the abstact, is obtained by taking all input parameters at their central 
values. It represents a large renormalization with respect to the initial value given by 
eq. (|2.5| ). In this respect, improving to 0(p A ) the chiral expansion, albeit challenging, 
might be needed in order to assess the degree of stability of the result. 



5 K L -K S MASS DIFFERENCE 

We apply the results of the previous section to the study of the Kl~K s mass dif- 
ference AMls- The full AMls can be split into short- and long-distance components 

as 

AM LS = AM SD + AM LD . (5.1) 

Notice that the "short-distance" component AM$d, generated by the lagrangian 
in eq. (|1 . 1|) , contains the hadronic parameter Bk- The value of Bk estimated in the 
previous part of the paper completes the determination of the box (or short- distance) 
component of AMls- 

We address now the issue of the evaluation of the genuine long-distance contribu- 
tion AMld- We will do it consistently with the evaluation of Bk- The interesting 
question is whether the interplay between AMsd and AMld reproduces the ob- 
served value AMr X q. 



5.1 Long distance AS = 1 induced contributions 



Many attempts have been made to estimate AMld [30] — 1 3T] by means of vari- 
ous techniques like chiral symmetry, dispersion relation with experimental data of 
s— wave 7T7T scattering, leading to a variety of numerical results. 

Our aim is to give a consistent estimate of AM L d based on the xQM approach. 
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NDR 
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GUQ2) 
G b LL (Qi) 
GUQ2) 

G 8(Q3) - 



(GG) 



G 8 _(Q,) = it 



GsiQs) = i <g> il ( 1 



M 2 

A 2 

x 



(gg) / 2 



M 



/? 1-6 



/\i 2 

A 2 

A 



Gll(Qi) 
GUQ2) 
G b LL (Qi) 
GUQ2) 

Gs(Q3) - 

Gs{Q±) - 
Gs(Q 5 ) - 
GsiQs) - 



~N~jt t 1 - 6 (GG) 

m 7T 
m 7T 

h-it i 1 ~ 6 (GG) 

\gg) 



'Ac: 



( 1 ~ 5 (GG) - 6^ 



it 1-6^ 



_2_ 1M 



/; 1 



9 f 



M 



/I 1-9 



iU 2 

A 2 

x 



Table 4: Values of the relevant AS = 1 weak chiral coefficients for two different 
regularization schemes: HV and NDR. The inclusion of the Wilson coefficients of 
the effective quark operators Qi is understood. 



As already mentioned, the K°-K° mass difference receives contributions from 
the exchange of the SU(3) meson field octet (we leave aside in this analysis the 
contribution of if) via the double insertion of the AS* = 1 chiral vertices. 

The complete bosonization of the AS = 1 lagrangian of eq. Ql.21 ) can be found 



in ref. P8| . Here we just quote the result for the operators Qi-e, which turn out 
to be relevant for the calculation of AM LD . The electroweak penguins Q7-10 gi ye 
a negligible contribution (of the order of 1%) due the smallness of their Wilson 
coefficients. 

The bosonization of the relevant operators leads to 

£aLi = G-8(Q3- 6 )Tr (\lD^D^) + 

G a LL (Q h2 ) Tr (AfEtZ^E) Tr (\l^D^) + 

G h LL {Q l>2 ) Tr (AfEt^E) Tr (\\Y?D»Y) , (5.2) 

where, as before, A* are combinations of Gell-Mann SU(3) matrices defined by 
(A})ife = 5u5jk and £ is defined in eq. Q3.2|). The covariant derivatives in eq. (|5.2|) 
are taken with respect to the external gauge fields whenever they are present. 

The notation for the chiral coefficients G^Q^q), G a LL {Qi^) and G b LL {Q\^) re- 
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a 



b 



Figure 5: One- loop long-distance contributions to the K° — K° mixing induced by 
the AS = 1 weak hamiltonian. The black box represents the insertion of the AS = 1 
chiral interactions. The octet mesons K, ir, rj are exchanged in the loop. 



minds us their chiral properties: G$ represents the (8 L x 1 R ) part of the interaction 
induced in QCD by the gluonic penguins, while the two terms proportional to G a LL 
and G\ L are admixtures of the (27 ' L x 1 R ) and the (8 L x 1 R ) part of the interaction, 
induced by left-handed current-current operators. These coefficient have been eval- 
uated in two different schemes of regularization HV and NDR, and the results are 
given in Table 4. In this table M is the constituent quark mass, that, consistently 
with previous analyses, we take at 220 MeV and A x is the chiral symmetry breaking 
scale (~ 1 GeV). 

The diagrams relevant to the evaluation of the long-distance contribution AM^d 
arise via one-particle and two-particle intermediate states (three-particle intermedi- 
ate states have been shown not to give significant contributions fl32"|). They contain 
two weak vertices, among those proportional to G§, G a LL and G\ L . Therefore we 
have to consider all the possible combinations: G a LL G a LL , G b LL G b LL , GgGs, G a LL G b LL) 

Using the Feynman rules reported in appendix A, it is found that the single 
particle intermediate state contribution give a result proportional to {Am 2 K — m 2 — 



3m^) which vanishes [30] by the Gell-Mann-Okubo relation. 

A non-vanishing contribution is obtained from the two particle intermediate 
states, which corresponds to the double insertion of the AS" = 1 chiral lagrangian 
as depicted in Fig. 5(a) and 5(b). To our knowledge, the relevance of the diagrams 
of the type (b) (tadpole diagrams) was first pointed out in ref. |$6|. 



The calculation is lengthy and the details can be found in refs. [40]. In eval- 
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A (4) 
ly QOD 


250 MeV 


350 MeV 


450 MeV 




0.113 


0.119 


0.125 


NDR 


Zl 


(0.0503) -0.524 


(0.0533) -0.663 


(0.0557) -0.781 




(0.982) 1.29 


(0.981) 1.39 


(0.980) 1.48 


z 3 


0.0180 


0.0360 


0.0870 


Z A 


-0.0471 


-0.0852 


-0.182 




0.0085 


0.0077 


-0.0129 


z e 


-0.0495 


-0.0947 


-0.226 


HV 


Zl 


(0.0320) -0.657 


(0.0339) -0.910 


(0.0355) -1.36 


z 2 


(0.988) 1.38 


(0.987) 1.58 


(0.987) 1.96 


z 3 


0.0137 


0.0301 


0.0798 


Zi 


-0.0292 


-0.0540 


-0.115 


z$ 


0.0070 


0.0100 


0.0123 


ZQ 


-0.0275 


-0.0515 


-0.112 



Table 5: NLO Wilson coefficients at /x = 0.8 GeV in the NDR and in the HV scheme 
(a = 1/128). The corresponding values at /i — m-w are given in parenthesis. In 
addition one has z 3 _ 6 (m c ) = 0. The coefficients Zi(ji) do not depend on m t . 

uating the loop integrals, we use dimensional regularization and modified minimal 
subtraction. 

5.2 AS = 1 Wilson coefficients 

In Table 5 we report the Wilson coefficients of the first six operators at the scale 
/i = 0.8 GeV in the NDR and HV 75-schemes, respectively. Since Re r in eq. ( |1.21| ) 
is of O(10~ 3 ), the K <-> K° transition is controlled by the coefficients Zi, which do 
not depend on m t . 
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5.3 AS = 2 Wilson coefficients 

The Wilson coefficients of the AS = 2 effective quark operator are denoted by 
7/i, 7/2 and 7/3 (see eq. QTTTD). 

The NLO calculations of 7/1 and 7/ 2 can be found in refs. and || respectively, 
while the analogous calculation for 7/3, which is particularly challenging, has been 
performed only recently by the authors of ref. [0. We have taken their results and 
evaluated the QCD factors for our choice of parameters. 

As an example, for Aqq D = 350 MeV, m[ pole ^ = 180 GeV, and /i = m c = 1.4 
GeV we find 

7/1 = 1.36 7/2 = 0.51 7/3 = 0.45 (5.3) 



5.4 Numerical Analysis 



According to Wolfenstein's notation ||39|| , we define the parameter D which char 
acterizes the long distance contribution 

AMld 

j-exi 
LS 



D = . (5.4) 

AMr 



Numerical estimates of D, AMld and AM so fc> r different values of Aq^, d and 
of the matching scale \x are given in Table 6. In this table we have fixed the gluon 
condensate to our central value (a s GG/7r) = (360 MeV) 4 and we have chosen for the 
quark condensate the value (qq) = —(280 MeV) 3 which gives us the best fit of the 
AI = 1/2 selection rule. The ranges thus obtained are 0.63 < AM L h s / AM^f < 1.11 
in the NDR and 0.54 < AMf s /AMf/ < 0.98 in the HV. 

A few comments are in order. Among the diagrams of Fig. 5(a), those containing 
two intermediate pion states dominate over kaon and eta exchange by about a factor 
of two. The diagrams of Fig. 5(b) (tadpole diagrams) give a contribution comparable 
in size with those of Fig. 5(a) but of the opposite sign, leading to a small and negative 
AMld in most of the parameter space. 

We disagree with ref. J5?J in the details of the calculation and on the the relevant 



interactions. In particular, the author of ref. [36] seems to neglect some of the leading 
insertions of the operator Q 2 . 

Our result depends sensitively on the value of the gluon condensate (the uncer- 
tainties in the short- distance coefficients related to varying the quark thresholds — in 
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Aqq D = 250 MeV 




li = 0.8 GeV 


fi = 1 GeV 




NDR 


HV 


NDR 


HV 


AM SD 


2.55 


2.42 


2.92 


2.80 


AM LD 


-0.34 


-0.38 


-0.35 


-0.38 


D 


-0.10 


-0.11 


-0.10 


-0.11 


AM th /AM exp 


0.63 


0.58 


0.73 


0.69 


A OCD = 350 MeV 




li = 0.8 GeV 


fi = 1 GeV 




NDR 


HV 


NDR 


HV 


AM SD 


3.07 


2.83 


3.56 


3.37 


AM LD 


-0.54 


-0.65 


-0.46 


-0.51 


D 


-0.15 


-0.18 


-0.13 


-0.14 


AM th /AM exp 


0.72 


0.62 


0.88 


0.81 


d = 450 MeV 




li = 0.8 GeV 


fi = 1 GeV 




NDR 


HV 


NDR 


HV 


AMsd 


3.91 


3.40 


4.55 


4.20 


AM LD 


-1.18 


-1.50 


-0.65 


-0.75 


D 


-0.34 


-0.43 


-0.18 


-0.21 


AM th /AM exp 


0.77 


0.54 


1.11 


0.98 



Table 6: Long-distance and short-distance box contributions to AMls, in units of 
10~ 15 GeV, for different values of the matching scale \i and Aq^ d in the xQM. 
We take for the gluon condensate the value (a s GG/7r) = (360 MeV) 4 and for the 
quark condensate (qq) = —(280 MeV) 3 , which are the values preferred by the fit of 
the AI = 1/2 selection rule at the same perturbative order. The "short-distance" 
component AM so is evaluated for a top quark pole mass of 180 GeV and for the 
values of B K given in Table 3. 
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Figure 6: The ratio (AM$d + AM^d) / AM^ P is shown as a function of the gluon 
condensate for Aq^ d = 350 MeV and fi = 0.8 GeV. We denote by (GG) the quantity 
(asGG/ir) 1 ^ in units of GeV. The dark and grey lines represent the HV and NDR 
results respectively. 



particular m c in the range 1.3 4- 1.5 GeV — affect AMls by less than 15% for our 
central value of Aqcd)- Fig. 6 shows the typical behavior for a choice of input 
parameters. The total theoretical mass difference AMjj S is a decreasing function of 
the value of the gluon condensate, analogously to the case of B K . If we let the value 
of the gluon condensate vary in the range of eq. ( |4.1| ), as we did in determining Bk, 
we obtain the overall range 

0.42 < AM%/AM™i < 1.40 (5.5) 

which represents our most conservative result. 

The scheme dependence of the result is satisfactory (< 20%) for most of the 
range of Aq^ d . Less satisfactory is the renormalization scale dependence. AMsd is 
an increasing function of \i and this behavior is not compensated by a corresponding 
decrease of the D parameter. This feature leads to a scale dependence that is about 
30% for Aqcd < 350 MeV - These results may indicate the need to extend the analysis 
to 0(p 4 ) in the chiral lagrangian expansion. An improved xQM calculation of the 



23 



AS = 1 and AS = 2 chiral coefficients at 0(p A ) is under way. 

A final comment about the width difference AT^s is necessary. A direct calcu- 
lation of the absorptive component of Fig. 5(a) gives about 1/6 of the experimental 
result. The reason is that the tree- level K — > tttt decay amplitudes do not reproduce 
the measured ones. Only by replacing in the vertices of Fig. 5(a) the one-loop 
results obtained in 0, we obtain the agreement with the experimental Ar^s. This 
is equivalent to computing directly the absorptive part of of Fig. 5(a) up to three 
loops. 

6 The Mixing Parameter lM\ t 

A range for the KM parameter Im X t , which is relevant for CP violating observ- 
ables, can be determined from the experimental value of £ as a function of B^, m t 
and the other relevant parameters involved in the theoretical estimate. 



Given m t , m c and the KM parameters 26 



\V US \ = 0.2205 ±0.0018 (6.1) 
\V cb \ = 0.041 ±0.003 (6.2) 
\V ub /V cb \ = 0.08 ±0.02, (6.3) 



we can solve the two equations 

(4) 
QCD 



s th (B K , \V cb \, \V US \, Agl D ,mi,m c ,77,p) = e exp (6.4) 



2 , 2 
T) +p 



w, 



12 



us 



Vr 



rb 



(6.5) 



to find the allowed values of the two parameters r\ and p appearing in the Wolfenstein 
parametrization of KM mixing matrix. 

We include in this analysis the interval of values for Aqq D in eq. (|1.9|) and for 
the gluon condensate the range in eq. (|4.1|) . The matching scale p is varied between 
0.8 and 1 GeV, while Bk is varied according to the range obtained in the previous 
sections. The results for m[ p ° l = 180 GeV are presented graphically in Fig. 7 . 

We can see that the equations ( |6.4| ) and ( |6.5| ) define two families of curves (re- 
spectively hyperbola and circles) in the (p-rj) plane. The allowed values of the two 
parameters correspond to the region delimited by the intersections between the two 
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Figure 7: Constraints on KM parameters from kaon physics. See the text for expla- 
nation. 



families of curves. The area enclosed by the two solid line hyperbolae corresponds 
to our most conservative range 0.54 < B K < 1.2. 

This procedure gives two possible ranges for 77 and consequently for lm\ t ~ 
^IKjsI |Kfe| 2 , which correspond to having the KM phase in the I or II quadrant (p 
positive or negative, respectively). For the central value of the top mass (jn t (mw) — 
183 GeV) we find 

< ImA, < 1.7 x 1(T 4 



0.67 x 10" 
in the first quadrant and 

0.41 x 10" 4 < ImA t < 1.7 x 10~ 4 



(6.6) 



(6.7) 



in the second quadrant. 

We also consider a "biased" estimate of Im At, obtained by fixing the gluon con- 
densate and Aq^ d to their central values and varying the matching scale fi between 
between 0.8 and 1 GeV. In this way we are spanning the following range for the 
renormalization group invariant parameter B^ : 



B K = 0.87 ±0.10, 



(6.8) 
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This choice of the input parameters restricts the hyperbolic band to the area 
enclosed by the dashed lines in Fig. 7. The overlap with the constraint of eq. ( |6.5|) 
leads to the following ranges of allowed values for Im At 



0.81 x 10~ 4 < ImA t < 1.6 x 10 



,-4 



(6.9) 



in the first quadrant and 



0.52 x 10 -4 < IniA t < 1.5 x 10 



,-4 



(6.10) 



in the second quadrant. 

Varying m t in the range m!f ole?> = 180±12 GeV affects very little the quoted lower 
bounds on Im X t while the upper bounds are changed by less than 20% (decreasing 
m t corresponds to increasing the upper limits). On the other hand, the upper bound 
on Im X t remains stable, beeing bounded by the maximum value of 77 obtained from 



We do not discuss here details of the bounds provided by AMg d which are repre- 
sented, for the central value of m t , by the area delimited by the grey lines in Fig. 7. 
The constraints of Bd — Bd mixing have a marginal impact in the determination of 
the overall range of rj. As the example in the figure shows, only the lower bound of 
Im X t in the second quadrant is affected by such an inclusion and raised to 0.6 x 10~ 4 . 
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eq. Q = 0). 
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A Feynman Rules for the AS = 1 Chiral Lagrangian 



We report the Feynman rules for the three terms (proportional to G a LL , G h LL and 
G£) of the AS = 1 chiral lagrangian which are relevant for the calculation of the 
long-distance contribution AMld- All momenta are entering the vertex. 



G a LL : 



K + {p l )^{p 2 ) 
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